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SECTION I. INTRODUCTION 


This report summarizes results from A study which had 
the objective of developing a vector wind gust model that is 
suitable for orbital flight test operations and trade studies. 
Detailed background information concerning earlier work can be 
found in References 1 and 2. In the work reported here, 
emphasis was given to verification of the hypothesis that gust 
component variables are gamma distributed, gust modulus is 
approximately Weibull distributed, and zonal and meridional 
gust components are bivariate gamma distributed. Section II 
describes a method of testing for bivariate gamma distributed 
variables; In Section III, two distributions for gust modulus 
are described, and the results of extensive hypothesis testing 
of one of the distributions are presented; Section IV establishes 
the validity of the gamma distribution for representation of 
gust component variables. Conclusions are presented in 
Section V. 


SECTION II. TESTING FOR BIVARIATE GAMMA 
DISTRIBUTED VARIABLES 


The hypothesis that absolute component gust and associ- 
ated gust length are bivariate gamma distributed can be tested 
according to the procedure described below. 

The probability density function of bivariate gamma dis- 
tributed variables is 


f (x,y;Yx'Yy=Y»P) = 


(3 3 

f(vra-p) 



EXP 



) 





2/p3j^3yXY 

P 


( 1 ) 


where I { } is the modified Bessel function of the first kind 
n 

of order n; 3 and 3„ are scale parameters; and y and y 
X y X y 

are shape parameters of the gamma distributions of x and y, 

respectively. The parameter y is the geometric mean of 

y and y . Estimation of y , y , $ , and 3„ from sample 
X y X y X y 

statistics is discussed by Adelfang (Ref. 1). 

Dimensionless variables T^ and T 2 are defined by 


= ^x^ 


T 


2 




( 2 ) 


The variables Tj^ and T 2 can be expressed in a coordi- 
nate system that is rotated by 45®; the transformed variables 
and Z 2 are given by 


2 



( 3 ) 


/I 

T 

(Tj+Tj) . 

/2 

T 

/7 

0 

(T,-T, ) • 

/T 

0 


(p^X + gyY) 
(pyY - p^X) . 


The probability density function of the transformed vari- 
ables is given by 





T-1 


1-p ^1 ^2 


(4) 


The probability that bivariate distributed variables 
and Z 2 will occur within the area bounded by the lines 
z^=zj, Zj^=Z 2 » and z^=-Z 2 (illustrated in Figure 1) can be 
calculated by numerical integration of the equation 




,/Jn e z^} dz^ 


= 


(1-p)^ (/2^) 


Y-^ 


r (Y) 


(5) 


where 


I„(W) = 


w 


n+2k 


k=0 


2nt2k p(n+k+i) 


n = Y - *5 


W 


_ /2p 

- r--^ 


( 6 ) 


Alternatively, P. can be estimated from the series; 


3 






A computer program for calculation of I , , using 

Equation (7) , has been developed. A sample of the calcula- 
tions of as a function of p for y = 3 is listed in Table 1. 
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OWGINAL PAGE IS 
OP POOR QUALITY 



standard results for checking the computer prograims are 
obtained from the closed-form solution for the case which 
can be expressed in the form 




1 



zf 1 


+ 1 ) 



1 

7p ‘ 


1 )] 


( 9 ) 


Values for are listed in Table 2 for selected values 

of zj and p . 


Table 2. P^(’f=l, p) Calculated from Equation (9) 



. . P 

z * 

. 1 

.25 

.50 

.75 

.875 

1 

.425980 

.445255 

.474470 

.495090 

.501805 

2 

.774586 

.774144 

.769772 

.763367 

.760084 

3 

.919308 

.911445 

.899446 

.889099 

.834464 

4 

.971975 

.965471 

.956084 

.948025 

.944361 

5 

.990370 

.986548 

.980820 

.975641 

.973206 

6 

.996704 

,994760 

.991624 

.988584 

.987097 

7 

.998874 

.997959 

.996342 

.994650 

.99'786 

8 

.999615 

.999205 

.998402 

.997493 

.99 /008 

9 

.999869 

.999690 

.999302 

.998825 

.998559 

10 

.999955 

.999879 

.999695 

.999449 

.999306 


Another useful special case is for p=0 and 2y equal to an 
integer. 


P^(p=0f 2y = an integer) 


= 1 


/2 

e 


Z * 

^1 


2y-1 

V 

k=0 


jn 


( 10 ) 
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The variation of P. as a function of correlation 

A 

coefficient, p, (for y=2) and as a function of shape parameter, 

Y, (forpsO.5) is illustrated in Figures 2 and 3, respectively. 

A comparison of observed and expected is illustrated 

in Figure 4; the line drawn at an angle of 45° to the abscissa 
represents perfect agreement between observed and expected 
values. Deviations of the plotted points from the line repre- 
sent differences between the observed and expected values. The 
data plotted in Figure 4 show a consistent pattern at 10 and 
12 km; for P^ <0.3, the observed is larger than the expected; 

for intermediate values (0.3 < P^ <0.8), the expected is larger 

than the observed. These results are the basis for initiating 
a more detailed amalysis of the validity of the gamma distribution 
hypothesis for the marginal distributiorts tcomponent gust and 
associated gust length) . The results of this analysis are 
described in the next section. 
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Figure 2. Series Approximation of P. as a Function of z* and p for 




































Figure 3. Series Approximation of P. as a Function of z? and y for 



■V‘ , 



A 



Figure 4 . Observed and Expected at 10 and 12 km 

Calculated from u Component Gust and Gust 
Length Data (A. = 2470 m) During February 

at Cape Kennedy 
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SECTION III. THE DISTRIBUTION OF GUST MODULUS 


Given that the absolute gust components are uncorrelated 
bivariate gamma distributed, then the probability distribution 
of gust modulus is obtained by numerical integration of the 
joint distribution expressed in polar coordinates. 


^lj2 


Pr{R < R*} = G(R*) = 


Y1+Y2-I 

nySrh.) ^ 

^ ^ 0 


[I] dR 


( 11 ) 


I 


Tf/2 “R(Bi cos0+B,sine) 

/ (COS0) (sine) ^ e ^ ^ de 

0 


where and B 2 scale parameters and and Y 2 the 

shape parameters of the u and v component gamma distributions, 
respectively. 

An expression that is approximately equivalent to 
Equation 11 is 

H(y. + Y,. AR*) 


where H(Yi + Y->» AR*) is the incomplete gamma function which can 

be calculated accurately with the series approximation given in 
Section II. 


12 



A 




l'*’^ 2 \ 

2r(Y^) r(>2) It-V*-) 


(13) 


Preliminary tests have indicated that reasonably accurate 
estimates of the probability distribution can be obtained from 
equation (12). However, it would be advantageous to determine 
if an alternative expression can be found which would not re- 
quire as much computation. The Weibull distribution, widely 
used in wind energy studies (Reference 3) was chosen to repre- 
sent gust modulus because of its relative mathematical simplicity 
and the availability of data for pvaramcter estimation. The cumu- 
lative probability function for the Weibull distribution of 
gust modulus is 


G(R*) = 1 - EXP 



(14) 


The paranK'-ters Ic 
approximation given by 


and c are calculated according to the 
Justus (Ref. 3) 


k 



--1.086 


(15) 


c 




(16) 


It is noted that equation (15) implies the relation. 


/'IbV = j^-1. 84162 , 
\ R / 
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whereas the exact relation for a Weibull distribution is 
given by 



r(l + 2 /k) 


T 

r^(i 


+ lA) J 


- 1 


( 18 ) 


The accuracy of the approximation has been evaluated for 
values of k from 0.5 to 10 by calculating the ratio, P, of 
the r: jht side of equation (18) to the right side of equation (17). 
Perfai t agreement is indicated when P=l. As illustrated in 
Tigure 5, for k > 1, P is within a few percent of unity; for 
k < 1, P approaches as k approaches 0. Therefore, it is con- 
ci uded that the approximation given by equation (15) is accurate 
for k '' 1 . The calculated values of k for gust modulus are 
between 2 and 3, which is within the range of acceptable accuracy 
of equation ( 15) • 


Parameters K and C, calculated from Equations (15) and (16), 
respectively, utilizing Cape Kennedy sample data are listed 
in Table 3. 

A comparison of the Weibull, the probability distribu- 
ticin associated with the modulus of a bivariate normal, and 
the observed • robability distribution is illustrated in 
Figure 6. It is indicated that there is little difference 
between the theoretical distributions for percentiles between 
20 and 98; for percentiles outside that range, the distribu- 
tions diverge; for this case, the observed distribution fits 
the Weibull slightly better than the bivariate gamma modulus 
distribution. 

The hypo basis that gust modulus at a reference altitude 
is drawn from a Weibull distributed population was tested for 
69 cases, 'i .e results summarized in Table 4 indicate that the 
hypothesis is accepted at the 0.05 level of significance in a 
large majority (65/69) of the cases. 
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Table 3. Pcuraraeters K and C of the Weibull Distribution for Gust Modulus at Cape Kennedy 
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Table 4. Summary of Results of Testing the Hypothesis* That 
Gust Modulus at a Reference Altitude (4, 6 ... 14 km) 

Is Drawn From a Weibull Distributed Population 
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SECTION IV. DISTRIBUTION OF GUST COMPONENT VARIABLES 


Four variables associated with gusts at a reference height, 
Hq, have been studied to establish the validity of the hypothesis 

that they are samples from gamma distributed populations. The 
four variables are illustrated in Figure 7. The variable is 

the largest u component excursion with sign equal to the sign of 
u at Hq; U 2 is the largest u component excursion of sign opposite 

Uj^ found by scanning upward after the second zero crossing associ- 
ated with Uj^. The vertical distance between Uj^ and U 2 is defined 
as L Range; the sum of the absolute values of u^^ and U 2 is defined 

as u Range. The variables u Range and L Range represent wind 
shear and wind shear altitude interval assc ^ted with gusts in 
the vicinity of Hq. Each of the four variables defined above have 

been calculated at six reference altitudes from a sample (150/ 
month) of February, April, and July Jimsphere wind profile data 
from Cape Kennedy. These data nets were tested to establish the 
validity of the hypothesis that each variable is drawn from a 
gamma distributed population. Acceptance or rejection of the 
hypothesis is at the 0.05 level of significance for a variate 
defined by 


'i 

X‘ 


N (0. 

T * 

i=l 



(19) 


= Observed frequency in the ith class interval 

= Expected frequency in the ith class interval 
(of the theoretical -gamma distribution) . 


The results of the hypothesis testing are described later. 


A. Absolute Gust Component and Associated Gust Length 


Gust, defined as the maximum excursion between successive 
zero crossings in the vicinity of a reference altitude, and 
associated gust length, defined as the distance between zero 
crossings, are each hypothesized to be drawn from a gamma* 
distributed population. The hypothesis is accepted at the 


1. Sample estimates for the parameters of the gamma distribu- 
tion are given in the Appendix. 
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0.05 level of significance, in a large majority of cases, for 
gust component magnitude (|u'|); specifically, the accept/ 
reject ratio is 47/22 and 46/23 for u and v con^onent magni- 
tudes, respectively. As indicated in Table 5, the ratio is 
significantly smaller for gust length (Lu and Lv) with rejec- 
tions exceeding acceptances (for method I) . The large number of 
rejections is attributed to large differences between observed 
and expected frequency of occurrence in the first few class 
intervals ; the observed frequencies are always much larger than 
the expected frequencies. Small gust magnitudes are associ- 
ated with small gust If.'gths that are observed as a consequence 
of the definition of gust used in this study. These small gust 
lengths are not measurable with the Jimsphere system; therefore, 
they are not considered to be valid data for hypothesis testing. 
By neglecting these data, we obtain the results summarized 
under II in Table 5 which indicate acceptance in a much larger 
proportion of the cases. 


B. U Range and L Range 


A summary of results of testing the hypothesis that the 
variables, u Range and L Range, are drawn from gamma o.’stributed 
populations is given in Table 6. It is indicated that the 
hypothesis for u Range is accepted at the 0.05 level of signifi- 
cance in 66 of the 72 cases. Acceptance is not a function of 
altitude except in July when the number of samples accepted at 
14 km was less than at the other altitudes. Acceptance was 
not related to filter choice with only slight exceptions 
(for A = 2470 during July and A = 6000 m during February 

w C 

one-third of the samples were rejected). Based on these results, 
it is concluded that u Range is gamma distributed. 

The results for L Range summarized in the lower half of 
Table 6 indicate acceptance of the hypothesis (46 of the 72 
cases) with not as strong a tendency as that indicated previously 
for u Range. Acceptance is an irregular function of altitude 
which is a minimum at 12 km where 50 percent is accepted to 
a maximum at 8 km where 75 percent is accepted. Acceptance 
is greater in July (75 percent) than in either April or 
February (58 percent for both months) . Acceptance is wecik or 
non-existent for A = 420 m and is strong for A large (2470 

and 6000 m) . 
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Table 5. Summary of Results of Testing the Hypothesis 
that u and v Component Absolute Gust and Gust Length are 
Drawn from Gamma Distributed Populations 


A/R 


( 2 ) 



Mathod 1 

nitar 

I 

11 . . 1 

X„(a) 

Month 

Month 1 


3 

4 

7 

All 

Non tha 

a 

4 
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Nontha 1 
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4/a 

6/0 

5/1 

15/3 

5/1 

6/0 


6/0 

17/ 

997 

6/0 

5/1 

4/a 

15/3 

6/0 

6/0 


5/1 

17 

3470 

5/1 

1/5 

a/4 

8/10 

6/U 

4/a 


5/1 

15, 

6000 

3/a 

5/0 

1/4 

9/6 

5/0 

5/0 


5/0 
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All Flltara 
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17/6 

ia /11 

47/33 

aa /1 

31/3 

a 1 /a 

64/5 


Lu 
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997 

3470 

6000 

5/1 

1/5 

0/6 

3/3 

4/2 

1/5 

1/5 

4/1 

4/a 

4/2 

3/3 

a/3 

13/5 

6/12 

4/14 

9/6 

6/0 

a/4 

4/a 

5/0 

5/1 

1/5 

4/a 
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16/3 
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32/37 

17/6 
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Flltar 

t 
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4 
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th« .OS laval of aiqnificarx* for varlata with a daoraai of fraadoai a - whara 

n • nuafaar of elaas tntarvals, b • nwbar of paraaatara of tha omm diatrlbutlon • 3. 

la tha ratio of tha nuabar of caaaa aceaptad to tha nwbar rajactad. 
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Table 6' Summary of Results of Testing the Hypotaesls 
That the Variables, u Range and L Range, at a 
Reference Altitude (4, 6, ... 14 km) are 
Drawn from Gaifiaa Distributed Populations 
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with ■ degrees of freedoD, a ■ n-l-b, where n - nunber of class intervals, b • pjadier 
of para*eters of the q a s ssa distribution ■ 2. 
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SECTION V. CONCLUSIONS 


This report has emphasized methods for establishing the 
validity of the hypothesis that observed gust varicJsles, 
including gust component magnitude, gust length, u Range, and 
L Range, have been drawn from gamma distributed populations 
and that observed gust modulus has been drawn from a bivariate 
gamma dist»*ibuted population that can be approximated with a 
Weibull distribution. An analytical procedure has been pro- 
posed for testi .g for the bivariate gamma distribution. The 
procedure has the advantage of not requiring frequency counts 
within narrow cells defined hy the intersection of intervals 
of tne marginal distribution; these frequency counts would be 
impractical and unreliable because of the limited sanple size 
(150) of the available data. Instead, the new method requires 
theoretical and observed frequency counting over larger areas 
associated with non-dimens ionali zed and transformed variables. 
Preliminary results utilizing this method have indicated larger 
observed than expected frequencies for small gust lengths and 
associated small gust magnitudes; this is attributable to the 
definition of gust used in this study. These small gust lengths 
are not measurable with the Jimsphere system and, as indicated 
in Section IV, the res\’lts of hypothesis testing for the margi- 
nal distributions are improved greatly by eliminating them from 
the date- sample. The hypothesis that gust component (u and v) 
magnitudes are drawn from a gamma distributed population is 
accepted at the 0.05 level of significance in 122 of the 136 
cases tested; for gust length (Lu and Lv) , 102 of the 136 cases 
were accepted. 

The variables u Range and L Range have been used to repre- 
sent component wind shear and shear interval associated with 
gusts. The hypothesis that u Range observations were drawn 
from a gamma distributed population was accepted at the 0,05 
level in 66 of the 72 cases tested; the acceptance ratio was 
somewhat smaller for L Range with acceptance in 46 of the 72 
cases tested. 

Testing of the hypothesis that gust modulus is drawn from 
a Weibull distributed population has yielded highly favorable 
results with acceptance of the hypothesis at the 0.05 level in 
65 of the 69 cases tested. 
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APPENDIX 


Parameters y and (J (calculated from sample moments) for 
hypothetical gamma distributions of gust component variables 
|u'|, jv'l, Lu, Lv, u Range, and L Range defined in Section IV 
are listed in Tables A-1 through A-6. 

The parameters in the tables can be used to derive 
the gamma probability density function of the form 


g(x) 


S£ 


Y-1 

rw ^ 


EXP(-Bx) 


(A-1) 


Equation (A-1) can be expressed in terms of a nondimensional 
variable y, i.e., y = , such that 


g(y) 


y^~^ EXP(-y) . 


(A-2) 


The probability that y does not exceed a specified value, 
Y, is given by 


y 

P {y < Y) - / g(y) dy 
^ 0 


1 ^ Y-1 

/ y^ EXP(-y)dy. (A-3) 


The integral on the right side of Equation (A-3) is the incom- 
plete gamma function, H(y,Y), which can be approximated with 
the series summation given by Equation 4 in Section II with 
the substitution 



a = Y 
X = Y . 
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Tad>le A-1. Gcunma Distribution Parameters y *nd £ of Absolute u Con^jonent Gust 
Estimated frora Sas^le Moment Statistics* 



6000 2.3696 1.5587 3.6039 1.8240 2.9507 1.4285 2.6261 | 1.1393 2.8570 .8762 3.0623 .7489 






Table A-3. Gamma Distribution Parameters y amd 6 of Aosolute v Component Gust 

Estimated from Sample Moment Statistics* 
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Teible A-4. rumma Distribution Parameters y and B of Gust Length, Lv 
Estimatftd from Sosqple Moment Statistics* 
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Table A-5. Geunma Distribution Parciraeters y emd g of u R£uige Estimated from Sample Moment Statistics* 
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Table A-6, Gamma Distribution Parameters y and 6 of L R£mge (Associated with U Range) Estimated 
from San;>le Moment Statistics* 
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